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When traditional linearised theory is used to study gravity-capillary waves produced by 
Q,' flow past an obstruction, the geometry of the object is assumed to be small in one or 

several of its dimensions. In order to preserve the nonlinear nature of the obstruction, 

asymptotic expansions in the low-Froude or low-Bond number limits can be derived, but 
Cd \ here, the solutions are waveless to every order. This is because the waves are in fact, 

exponentially small, and thus beyond- all- orders of regular asymptotics; their formation 

is a consequence of the divergence of the asymptotic series and the associated Stokes 

Phenomenon. 

^ ^ In Part 1, we showed how exponential asymptotics could be used to study the problem 

^\ • when the size of the obstruction is first linearised. In this paper, we extend the analysis 

Cn \ to the nonlinear problem, thus allowing the full geometry to be considered at leading 

^^ ■ order. When applied to the classic problem of flow over a step, our analysis reveals the 

existence of six classes of gravity-capillary waves, from which two share a connection 
'nT . with the usual linearised solutions first discovered by Lord Rayleigh. The new solutions 

^~^ \ arise due to the availability of multiple singularities in the geometry, coupled with the 

interplay of gravitational and cohesive effects. 

Key Words: surface gravity waves, capillary waves, wave-structure interactions 

1. Introduction 

Consider water that flows pa st a fishing line , whic h we model as a pressure distribution 
applied near the surface. As iLord Ravleighl ( 18831 ) demonstrated, if the velocity of the 



stream is kept above some critical velocity, then capillary waves are produced upstream 
and gravity waves downstream. However, if the speed of the stream is too small, then 
subcritical or supercitical solitary waves are produced instead. This classification of the 
dynamics holds similarly for many flows past more general obstructions, where we can 
linearise the geometry in one or several of its dimensions. Thus, typical linearised theory 
does not really distinguish between flows pas t differently s hape d objects, but mere ly 



requires that they are sufficiently small [see e.q. lPorbed (|l983l ) and lKing fc Blooij (|l987l) ]. 



The question which we address in this paper is: what can he said about gravity- capillary 
flows over obstructions which are not linearised? 

There are two important parameters: the Proude number, F (ratio of inertial to gravita- 
tional forces) and the Bond number, B (ratio of gravitational to surface tension forces) . In 
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a previous work (jTrinh fc Chapmadl20136l — henceforth referred to as Part 1) we showed 
how the typical hnearised theory can be re-interpreted in the fow-Froude and low-Bond 
limits. If q is the speed of the free-surface flow, and F^ — 0{e) while B — 0{£^), then 
the typical asymptotic expansion gives the base series, which we denote as ®, and write 



® 



qo + egi + e^92 



Oie') 



(1.1) 



valid as e — )■ 0. The base series, ®, is waveless to every order, and in fact, the gravity 
and capillary waves are exponentially small and thus beyond-all-orders. These waves are 
switched-on when ® is analytically continued across critical curves {Stokes lines) in the 
complex plane in a process known as the Stokes Phenomenon. We use the notation ® > 
© or ® > © to indicate the switching-on of a gravity or capillary wave, respectively, 
and these waves are written as 



Ai+0{e) 



= -X9/e 



and 



© 



Bi+0{e) 



-XcA 



(1.2) 



In the linearised problem of flow over a step (Part 1), Stokes lines originate from the 
single, merged singularity which represents the step, and the point where the gravity 
and capillary waves are equal in magnitude corresponds to the critical bifurcation in the 
Froude-Bond plane. 

The nonlinear analysis of this work differs from the linear analysis in two important 
ways. First, a nonlinear geometry usually contains multiple singularities and hence mul- 
tiple Stokes lines; each singularity then has the potential to produce both gravity and 
capillary waves. Second, gravity and capillary waves can themselves interact, with © > 
© or © > ©. These more complicated secondary switchings may be accompanied by 
crossing Stokes lines, the higher-order Stokes Phen omenon, and other, more advanced, 
aspect s of exponential asymp t otics highli ghted by Howls e.t al. (|2004l ). lOlde DaalhuisI 
(|2004l ).IChapman fc Mortimeij (|2005[) , and lChapman et al\ (|2007h . 

Our main result shows that for a nonlinear obstruction, the typical bifurcation curve in 
the Froude-Bond plane thickens, revealing new possibilities for solutions. The thickness 
of this curve is a manifestation of the finite nature of the obstruction, and shrinks to zero 
as the geometry is linearised. For exa mple, in the case of a rectangular step, rather than 
the two standard linear solutions of iLord Ravleighl (|l883^ . there are now six possible 
solutions, for which the previous two are special cases. 



1.1. Physical and mathematical background 

A review of the physical and mathe matical literature that mot ivates this work can be 
found in the introduction of Part 1 (|Trinh fc ChapmanI 1201361 ). The situation of free- 
surface gravity-capillary waves produced by flow past an obstruction is a well studied 
problem, and so we s hall only mention the more comprehensive reviews of the to pic by 



^ , ^ p 

Dias fc Kharil (119991) an d IVanden-Broeckl (|2010[ ). and of course, classic texts by iLamb 
(|l932h and lStokeil(ll957l) . 



As for the mathematical techniques we use in this paper know n as exponential a sym.p - 



totics, these are directly based on the methods introduced by IChapman et all (|l998|) 



and previously applied to study wave-structure interactio ns for pure capillary waves 
(jChapman &: Vandcn-Broeck 2002), and p ure gravity waves ( Chapman fc Vanden-Broeck 



20061 iTrinh et a/.ll201l[ iLustri et a/.ll2012l ). A variety of other app roaches to exp o nential 



asymp totics are available, and we refer the readers to the books bv lDinglg (|l973[ ). ICostin 



asymp totics are avanaoie, ana we reter tne reaaers to tne dooks pv iumgiei ( 
(|2008[ ). and iBovd (, 19981 ) for an overview of the of the different techniques. 



D 
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Figure 1. We consider a typical flow over an obstruction. The flow in the physical xy-pla,ne 
(left) is flrst mapped to the potential w = (j> -\- itp plane, then again mapped to the upper-half 
C, = £, +irj plane (right). For flow in a channel, the latter map is given by (^ = e~™. 



2. Mathematical formulation 

We briefly recapitulate the relevant equations, which parallel the ones presented in Part 1, 
but this time, we allow for the possibility of flows over more general geometries. Consider 
steady, two-dimensional potential flow of an incompressible fluid with upstream velocity 
C/, and a prescribed length scale L. The flow is non-dimensionalised with characteristic 
scales of U and L/tt for the velocity and lengths, respectively, and the physical z — x + iy 
plane is mapped to the complex potential w — 4> -\- iip plane, then mapped again to the 
upper- half C, = S,-\- irj plane using Q — e^™. The free-surface is then given by 



\ogq 



de, 



/3e 



9^$' 



l3Te^ 



1 r o{e) 

2(f9 dqdO' 
1 "777 +9-77-77 



(2.1a) 
(2.1&) 



where F^ — f3e is the square of the Froude number, B — /3re^ is the Bond number, and 
for hodograph variables \og{d'w / dz) — log q — i9, where q is the fluid speed and 9 is the 
angle the streamlines make with the a;-axis. This set-up is presented in Figure [T] 

Our analysis must then be extended to the complexification of the free-surface (orig- 
inally, q and for ^ € M+). Complexifying the free-surface into the upper- half ^-plane 
and relabeling ^ H- C and cj) ^-^ w gives 



/3e 



logq- 
^dq 



9{a 



dw 



Pre^ 



00 s s 

d'^9 dq d9 



de+.^9{c), 



dw'^ 



q 



dw dw 



sm9. 



(2.2a) 
(2.2&) 



where J^ denotes the Hilbert transform operator on the semi-infinite interval ^ ^ 
which corresponds to the free-boundary, i.e. 



.Jif9{C) 



1 



mi 

e-c 



de. 



(2.3) 



2.1. The inclined step 

The methodology developed throughout this paper is applicable to most general non- 
surface-piercing geometries, which are chosen by imposing the value of 9 along the neg- 
ative C = ^ axis. With minimal modification, it can also be applied to surface-piercing 
obstructions (such as for a ship) , th e only signifi c ant d ifficulty being the assumptions-to- 
make near the point of contact [see lTrinh et al\ ( 20111 )]. 
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However, for concreteness and illustration, we will often take as an example the step 
in a channel, with 

C = e-™ (2.4) 

and 

for C < -b 

forCe(-6,-a) (2.5) 

for C e (-a, 0) 

where < a < 5. We will always imagine a stepping up (from Icft-to-right), with b 
associated with the stagnation point, a associated with the corner, and therefore < 
fT < 1. In the next section, we will see that at zero Froude and Bond numbers, the 
leading-order solution is provided by substituting the above values for 9 along ^ < 0, as 
well as = for C > 0, into the boundary integral equation (I2.2flp . giving 

4^1 and 61 = + 0(6), (2.6) 

C + aJ 

which is often termed the rigid-wall flow, as it is equivalent to replacing the free surface 
by the rigid wall, 6* = 0. 



3. Asymptotic approximation 

Substituting the usual perturbation expansions (the base series), 

c» oo 

0--^e"^„ and -? ^ ^ £>„ (3.1) 

into (p:^ and ((236|) yields at 0(1), 

00 = 0, (3.2a) 

loggo = -- r ^de, (3.26) 

and at 0{e), 

go^ = -/i, (3.2c) 

^-iei = j^di{c). (3.2d) 

Notice that imposing the bottom topography in (J3.2&P determines 9 along ^ e M^. The 
leading order solutions (|3.2ap - (j3.2&p . simply correspond to the rigid wall flow whereby 
the free-surface is replaced by 6* = 0. 

The full expressions for the higher 0{e") terms are prohibitively complicated, but since 
we are mainly concerned with the limit n — )■ oo, we may proceed in the following manner: 
the complexification of the leading-order free surface, go will typically contain singular- 
ities, identifiable with singularities in the flow-domain, such as those corresponding to 
corners or stagnation points. However, because all the higher-order problems are linear, 
no new singularities can be introduced and thus, for all n, the singular points of (?„ must 
be those same points as for qq. 

Now if we examine the dynamic condition (J2.2&P , we can see that each successive term 
of the asymptotic approximation requires the derivative of the previous term. Then if (/„ 
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contains a singularity of the form 1 / (w — w*)^ , Qn+i will contain a singularity of the form 
n/{w — w*)"+^. Thus as n — >■ oo, we can expect the late terms to behave like factorial 
over power, or 



e(w)r(n + 7) 



and 



g« 



g(w)r(n + 7) 



(3.3) 



In fact, the inductive argument used to explain the divergence of the asymptotic ex- 
pansion according to (|3.3p also serves to explain why the free surface waves are expected 
to be exponentially small in the limit e — ^ 0: the leading order solution (j2.6|) is waveless, 
and since the values of g„ and On only depend on the derivatives of the previous orders, 
then none of the terms in the naive expansion (I3.ip will contain waves. 

Proceeding then with the ansatz p.3p . we can now pinpoint the necessary terms re- 
quired at 0(e"). In the limit that n — ^ cx), terms like qmin for m finite dominate terms 
with smaller indices in n, such as qmQn-i- Moreover, differentiating a term increases the 
order (in n) by 1, so a term like edqn-i/dw is of the same order as (?„. The relevant terms 
at 0{e") of the boundary integral equation (|2.2a[) are thus 

Q7I (I7I-IQI -Q ^^Q r^s 

90 9o 

for n ^ 2. It is known that as n — > cxd, the Hilbert transform on the right-hand side of 
this equation is exponentially subd ominant to the terms on the left [see Part 1 or e.g. 
Chapman fc Vanden-Broeckl ( 2006f ) for further details]. Thus, 



<7o 



iqiqn-i 



00, 



(3.4) 



or. 



Qn ^ iqo^n + ^^n-l?! + . • . aS TT, — !■ OO. (3-5) 

As for Bernoulli's Equation (j2.26p . we will use p.Sp to replace g„ with 0n, after which 
the relevant terms at 0{e") are 



first and second order as n — ^ oo 



^^O 


O'n-l - 


^9o 


01 


-2 + 


1 


On 








+ 


Siq^c 


-I'o 


On-1 + 



iiqlqi - Tqoq'g 



9' 

^n-2 



2Tqoqi 



:-3 = 0, (3.' 



second order as n — > oo 



with primes for differentiation in w. Using the ansatz p.3|) in the above equation, we 
have at leading order as n — >■ oo. 



which is simply solved to give 



^q'oX'-rqUxr + ^ = 
dw '^Tqo 



where 



4t 



A = q^- A and ^ = ^ 

P 



(3.7) 



(3.8) 



(3.9) 



and A will turn out to be a key parameter. Remember that x is the portion of the ansatz 
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(j3.3p that expresses the singularities of the higher order terms. Since x(u'*) = for one 
of these singularities, w* , we may express 



X±iw) 



9o±Va' 



Srqo 



dip. 



(3.10) 



Note that in writing (|3.10|) , we shall restrict the path of integration to be along the same 
Riemann sheet as w* . Integration contours that traverse different sheets (for example, 
by crossing the branch cut(s) from A = 0) will be addressed in later sections, and in 
particular. Appendix [Xj By t aking the lim it of (|3.10p as r — )■ and comparing with 
Chapman and Vanden-Broeck ( 2002t 120061 ). we see that the positive sign corresponds to 
capillary waves and the negative sign to gravity waves. 

Now returning to the dynamic condition (j3.6p and with the ansatzes (J3.3I) . we find at 
next order in n: 



ill 



& 



^lo 



[-2x'e'~x"Q}+[^iqWo 
iiqlqi - TqoQQ 



We may write this as 



with 



+ 



e 



e 



1 A^ A' 

"4 a" 2y^A{A + AJ 



2Tqoqi 



F±iw), 



{(x')'e} = o. (3.11) 



F±{w) = - 



1± 



A + 2A 
v/A(A + A) 



and thus, after one integration, 

(^/ATA + ^/A)^'A± 



e(w) = 



Ai/4 



X exp 



F±i(p) dip 



(3.12) 



(3.13) 



(3.14) 



where A± is a constant and we may begin the integration at any arbitrary point w = w* 
where the integral is defined (often w* is a natural choice for w* , but the integral may 
not exist at w*). Note that p. 51) . allows us to relate Q and Q, with 



Q ^ iqoQ. 



(3.15) 



As a check, we may take the limit of r — > and use the negative sign of (|3.14p . This 
recovers the pre-factor for the gravity-wave problem. 



e_(w) 



A_i 



'3i 



— dip 



* % 



(3.16) 



to leading order, which differs from (3.16) in IChapman fc Vanden-Broeckl (|2006f ) by a 
factor of a half (our A_ is their 2A). 

A note on the ± notation 

As we have seen, each singularity is associated with two branches of x = X± in (I3.10|) . for 
which the positive sign has been shown to correspond to capillary waves and the negative 
sign to gravity waves. There are also associated 6 = 0±, Q = Q±, and A = A± quantities. 
In what follows, we will sometimes use non-subscripted variables if it is unimportant 
which branch of x we are referring to. 
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4. Optimal truncation and Stokes line smoothing 

The underlying divergence of the asymptotic expansions will cause the Stokes Phe- 
nomenon to occur: as the complexified asymptotic solution crosses a critical line (the 
Stokes line), a small exponential switches on. Because the switching-on of the exponen- 
tial is almost always via an error function ( Berrvlll989f) . the optimal truncation and 
Stokes smoothing procedure will be similar to the one we performed in Part 1, with the 
exception of the nonlinear terms. 

To begin, we truncate the asymptotic series at n = A/" so that 



AA-l 
n=0 



A/'-l 



Rat 



and 



q = 



n=0 



'q,i + S^f, 



where the remainders are related by p.Sp and thus 
Sat eqiSjsf 



e^gi^A^- 



90 



ll 



ql 



iRj\f. 



(4.1) 



(4.2) 



We will substitute the truncated sums (|4.ip into Bernoulli's Equation (J2.2&I) and in doing 
so, we will see two separate types of terms, the ones involving only g„ and 0„, and the 
ones involving the remainders, Rjsf and Sj^. 

Let us first study the terms involving the remainders. The remainder 5V can be written 
in terms of Rj\f by (|4.2p and after making this substitution, we are particularly interested 
in the leading-order terms, which are indicated by factors of Rj\f, eR'j^, and e^R'^, and 
second-order terms, which are indicated by factors of eR^/, e^R'j^, and e^R'J^. The relevant 
terms from Bernoulli's Equation are then given by the linear operator defined by 



C{R^; e) 



^Qo 


6i?V- 


tQq 


e'R'^ + 


1- 


Ru 




+ 


3iq^q'„ 


eRj^ + 


Siq'^qi 



TqaQo 



2 jyl 



'A/' 



2t(7o<7i 



e^i?^, (4.3) 



which is exactly the same form as the left of p.6p . We then introduce the Stokes smooth- 
ing parameter, S — 5(w), and set R_\f = 5[0e~^/^]. We know that the ansatz Qe^^^'^ 
solves £ = 0, so only terms involving derivatives of 5 will be left. After some computation, 
we find to leading order. 



CiR^;e)^eQe-^/' 



dS_ 
dw 



ill 



2rqlx' 



Writing dS/dw = x'dS/dx, and using (|3.7p gives 
C{R^f■, e) 



eOe-^/'-f. 
dx 



-,+rqlix'f 



Now let us turn to the terms involving g„ and On'- when the truncated sum in ()4.1|) is 
substituted into (j2.2b| . terms of ©(e ) will automatically be satisfied, and this process 
leaves us only with the remnant ©(e-'*^) contributions from the inertial terms, as well as 
the 0{t ) and C(e-^+^) contributions from the surface-tension terms: 



e^ 



2dqjsf-i 



dw 



e^r 



2 d 6'a/-_2 

'^" dw^ 






e-^+V 



2d ^AA-l 



% 



dnp' 



e^+V 



diiP- 



(4.4) 
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and thus in total we have 



C{R^f:e)r^e 



A/" 






2d 9_\f-i 



dw'^ 



(4.5) 



Since we are interested in the hmit e — > when A/" — > 00, we can substitute the late-orders 
ansatz of p.3p into the right-hand side of (|4.5p . giving 



C{R^; e) 



X 



■J^+i 



l-^rAx'?'-^^^^ 



re 



ii) 



(4.6) 



The optimal 



We introduce a coordinate system along the Stokes line using x 

truncation point (where adjacent terms of the expansion are equal in size) is at [r/e] , so 

we write A/" = r / e + p where p e [0, 1). Changing to differentiation in i9, using 



and applying Stirling's formula to 



dS_ 

d-d 



rqlix'f 



d 


.e-'^ d 


" r di?' 


mula to (|4.fi|) gives 


\/27rr 

£7 + 1/2 


[1 , r^A^?\ 

_/3 e'« 



_j^Njr/e+p+7 ^re'"'/€g-r/e 



(e--) 



The exponential factor on the right is exponentially small, except near the Stokes line 
"(5 = 0, where the critical scaling occurs with ■& = \fl'Q. In total, then. 



dS 



'27rr 



e^ 



■ exp 



r-d^^ 



and we have recovered the typical error function expression for the Stokes multiplier, S. 
We now integrate this expression across the Stokes line from upstream, or i? = 00, to 
downstream, or i? = —00. The apparent jump in the remainders of 9 and q are then 



Rn 



S^f 



downstream 

upstream 
downstream 



27ri 



-ee 



-xh 



£7 ^ 



^cxp 



(4.7a) 
(4.76) 



upstream 

We can re-introduce the notation for the switching-on mechanism used in Part 1. We 
write, for example. 

Stag. 



E^>" 



-^ Ve"9n 



Qcxp 



^ ^ ^" ® >© 



Corn. V~^ n 

> 2^ e qn + gcxp, 



and the arrow notation should be read as "the base series turns on a capillary /gravity 
wave as the Stokes line from the stagnation/corner point is crossed" . 

To finalise the analysis, we need to also complexify the free boundary into the lower- 
half plane. This analogous process yields the functional complex conjugates of (|4.7aP and 
(|4.76p . and thus the total leading order contribution along the free-surface is given by 
twice the real parts of (I4.7ap and (|4.76p . 

gexp, total - -^S (Qe^^/^) , (4.8a) 

^exp,totai^-^3(ee-x/^). (4.8^,) 
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In order to fully determine the waves, we must also compute the values of the pre- 
f actor, A, that appears in p.l5|) and p.l6p for Q and 9, and also the value of 7 in ()4.8p . 
This can be done by matching the outer solution (I3.ip with an inner solution valid near 
each of the singularities. Although this was easily done for the linear p roblem in Pa rt 1, 
the nonlinear problem is much more tedious, and readers can refer to iTrinhl ( 20101) for 
the detailed procedure. 



5. Inner limits of x and Stokes lines 

In the low-Froude and low-Bond limits, we can think of the exponentially small free- 
surface waves as having been generated by singularities in the flow domain. For example, 
in the case of the step (|2.6p . both the corner and stagnation points are respons ible for 
produ cing waves through their associated Stokes lines. Furthermore, we know from lDingld 
( 19731 ) that Stokes lines with ® > © or (g) > ©, are given at points where the base 
series (with phase zero) reaches peak exponential dominance over the capillary or gravity 
waves (with phase — x/e). This yields the two conditions 



%X±) = and 3?(x±) ^ 0, 



(5.1) 



and for a given geometry, the trajectory of the Stokes lines can be derived by numeri- 
cally evaluating the integral p. 101) and applying conditions (|5.1I) . Whether a Stokes line 
encounters the free surface is a problem that thus depends on the global behaviour of q^; 
however, we can still study the local properties of the Stokes lines near their associated 
singularities. 

5.1. Local Stokes line analysis for general channel flows 

Let us assume that there is a singularity in the outer solutions p.ip at the points z — z* , 
w — w* , and C = C*i which correspond to the different planes introduced in Figure 
[1] We introduce the shifted coordinate W = w — w* , and then from potential theory, 
W ^ const X (z — z*)'^ for some k; in the case of stagnation points, k = 2, while for 
corner singularities, k — t:/v, where v is the in-fluid angle of the corner. In the analytic 
continuation of the free boundary, we have from the definition of the complex velocity, 

^'^ ^-'^" =qo = const X (z - z*)"-'^ = const x Ty(«-i)/«, 



dz 



Qoe 



(5.2) 



and thus the inner limit of the outer solution is given by 

g-c(u;-u;*)" = cW°' 



(5.3) 



where a — {k ~ 1)/'* and c is constant. 

Using X in (13. 8p and the limiting form for qq in (15. 3p . we see that x must exhibit 
different limiting behaviours as the singularity is approached, depending on which sign 
of the square root is chosen and whether a is positive or negative: 



X± 



ic 



i 
' c^l3{l-Za) 



T- 



le 



Arg(A)/2 






ATg VF"+i for capihary {+) and a < 



X^W 



x^w^-" 



1 3" for gravity (-) and a <0 (5.4) 
for a > 



cv^(l - a) 
For a < 0, the choice of the negative sign for the square root leads to gravity waves. 
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whereas the positive sign is associated with capillary waves; this is clear from the ap- 
pearance of either /3 (gravity) or r (capillary) in the coefhcients of Xi but can also be 
understood by taking r — ?► 0, and noticing that the capillary root in p.7|) disappears to 
infinity, leaving only the gravity root. 

There is a notable difficulty in determining the local behaviour of x± in (|5.4p when 
a > 0, and this deals with the value of Arg(A), which must tend to tt or — tt as we 
approach the singularity from the upper half-C-plane. For the purpose of the Stokes line 
analysis in the next section, it is sufficient to see that if we approach the singularity from 
along the upstream boundary, and set W — \W\e^^ , then ^{qq — A) tends to zero from 
below as \W\ -^ 0, and thus 

gArg(A)/2 ^ _ .^ (5 5) 

We also mention that the assumption that the potential follows a power of {z — z*) in 
(|5.2p will not be true for the case of a sink or source, where W ^ const, x log(z — z*). 
However, this analysis requires a somewhat different approach because the location of the 
singularity corresponds to \w\ — ^ oo. The problem of applyi ng exponent i al asy mptotics 
to gravity flow past a submerged line source was studied bv lLustri et al\ ( 20131 ). 



5.2. Local Stokes line analysis for the step 

The existence of a Stokes line does not necessarily guarantee the appearance of surface 
waves. Indeed there are cases where the associated Stokes line fails to intersect the free- 
surface, as well as cases where the line lies on a different Riemann sheet altogether. 
In this section, we study the simplest scenario of a Stokes line which emerges from 
its singularity in the 'in-fluid' direction of the upper-half ^-plane. Whether such a line 
actually encounters the free-surface (C G R"*") is a global function of the leading-order flow 
and in general, must be checked by evaluating the integral in (|3.10p . For the moment, 
we ignore lines which leave in a direction of 'into' the boundary (3(C) < 0), and lines 
which enter secondary Riemann sheets. From a study of the global contours of x, we have 
observed that relevant interactions are generally governed by Stokes lines that traverse 
simple paths (from singularity to free-surface along the same Riemann sheet). However, 
in ^ and Appendix \^ we shall see the role played by more complicated Stokes line 
trajectories which cross onto different sheets. 

Recall that the value of a in (|5.3p determines the boundary's geometry near the singu- 
larity, with a < for a corner and a > for a stagnation point. The value of c determines 
the orientation of the geometry, and this is not quite arbitrary because gravity provides 
a reference direction. Again, let us consider the two singularities associated with the step 
geometry in ()2.6|) . 

For the case of the corner, with a E ( — 1,0) and c e R+, then in the notation of ()5.4p . 

J |X3|e^"/2^"+i for capillary {+) 

^^^ ^ |X2|e-"/2^i-3" for gravity (-) 

as W — )■ 0. If we write z9grav = Arg(W) for the local angle of a gravity Stokes line, and 
^cap for the local angle of a capillary Stokes line, then we have from (|5.1I) that 

2m+l/2\ ^ „ f 2m +1/2 

— — and ?9cap = vr — — 

l-3a/ \a-fl 

where m g Z. As a increases from —1 to 0, the first (to — 0) gravity Stokes line increases 
from i?g,av = tt/S to 7r/2. However, for the same range in a, the first capillary Stokes line 
only enters the upper-half plane beginning at a = —1/2, where ^cap = tt. As a increases 
to 0, t^cap increases to tt/2. By taking the local analysis to next order in W, it can be 
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shown in the case of the Stokes hne initially leaving along the boundary, i?cap = vr, it 
must continue along the boundary until it encounters the stagnation point. 

For the case of the stagnation point, with a E (0, 1) and c = |c|e^'^'", then (|5.4p gives 

X±-T\Xi\e^'''W^-", 

as W^ — >■ once we have used (|5.5p . The Stokes lines leave at angles of 

^Eiav = TI" ( — 1 and l^caD = TT 



grav - M ^_^ J c..^ -cap - ^ ^ _ ^ 

for TO G Z. For gravity waves, there are no relevant Stokes lines lying in the upper half- 
plane. For capillary waves, the only relevant value of tti is 1, which shows that there 
exists a Stokes line along 6 = tt for all shapes. However, a secondary analysis to next 
order in W shows that the Stokes line which initially tends along ^cap — tt does indeed 
leave the boundary into the upper half-plane. 

The local Stokes line angles in the physical plane can be retrieved by using (|5.2p and 
(|5.3p . If z = z* is the corresponding singularity, then 



1 



c(l-a) 



W^-". (5.7) 



In summary, given a step-up obstruction with initial inclination cttt, the following 
occurs: first, the stagnation point always produces a capillary Stokes line; second, the 
corner always produces a gravity Stokes line, and only produces a capillary Stokes line 
if (T > 1/2. Table [1] provides an illustration of our local Stokes line analysis. Note that at 
a = 0, there is no singularity, so the corresponding entry instead refers to the limiting 
process of a ^ 0~. 



6. Turning points and Stokes lines for the step 

The singularities encountered in ^are not the only critical points in the gravity-capillary 
problem. From p.9p . the turning points, C = Ct, found at 

A^q^-A^q^O-J = 0, (6.1) 

also represent a type of wave-generating singularity. These points do not share a connec- 
tion with the physical geometry, but rather, they correspond to locations where deriva- 
tives of the capilla ry and gravi ty wavenumbers are equal. By drawing an analogy to 
the Airy equation (jWhitd l2005l Chap. 11), we would expect that these turning points 



can also produce Stokes lines, across which gravity and capillary waves can themselves 
interact. 

Care must be taken in solving A = 0, however, since this generally involves the inver- 
sion of composite complex- valued functions with multiple branch cuts. In the case of the 
step flow (|2.6p . for example, since A > 0, all the possible solutions are given by 

^^ + M =^i/4>.e^"''/^". (6.2) 

Cr + aJ 

Different values of A: G Z in (j6.2p may correspond to different turning points, but we must 
always check (a posteriori) that the final location lies on the main Riemann sheet {i.e. 
where qo is real and positive on C G R"*"). Like the previous Stokes line analysis of §5.2|, 
we only consider points on this immediate sheet; numerical computations of the global 
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Local power (a) Grav. Stokes (i9grav) Cap. Stokes (i9cap) Physical sketch 
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Table 1 . A summary of the local behaviour of a Stokes line (shown dashed) as it emerges from 
a singularity for which the complex velocity behaves like dw/dz ~ c{w — w*)" . Angles & are 
relative to the positive 5R(i(;)-axis. Angles in the physical plane follow from H5.7I) and the solid 
arrow indicates the direction of flow. 



contours of x seem to indicate that critical points that lie on secondary sheets do not 
affect the free surface (though there are exceptions: see Appendix R)) . 

There are three important turning points in (|6.2|) . which we denote as Ct = Cii C2, and 
^3, corresponding to fc = 0, 1, and —1, respectively. Their locations for various values of 
a and A are shown in the main portion of Figure [5J Recall that the upper half-C-plane 
shares a connection with the z-plane, and we have illustrated the physical connection 
within the insets (a) to (f). Note that we have taken the branch cut, marked by a striped 
line, along the real axis in the ^-plane, between the corner and stagnation points. 

First consider the two complex turning points, C2 and Cs- For a fixed values of A and 
(T > 0.5, the point C2 lies somewhere above the real axis, and Ca is its conjugate. For fixed 
a and as A — )■ 0, C2 tends to the stagnation point, while \i A ^^ cx), C2 tends to the corner 
point. This motion for variable A is shown as a solid curve in Figure [2] On the other 
hand, if A is fixed and cr — > 0.5, C2 moves towards the branch cut, so that at cr = 0.5, ^2 
lies directly on the angled step, while C3 has moved onto the adjacent Riemann sheet. If 
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Figure 2. Movement of the turning points, C,\, C2, Cs lor the step with 6 = 2 (stagnation point) 
and a = 1 (corner). Exact values are shown in the ("-plane and the insets provide depictions in 
the physical plane. The C,\ points are black squares, and (2,3 are grey squares. For fixed values 
of A and decreasing values of a (thus fixing the flow parameters, but changing the geometry), 
the different values of (2 and (3 are shown as a dashed line in the main illustration and through 
insets (a) to (c). Similarly, for changing values of A and fixed values of a (thus changing the flow 
parameters, but fixing the geometry), the values of (2,3 are shown as a solid line in the main 
illustration, and through insets (d)-(f). Note that for a < 0.5, both (2 and ("3 have passed onto 
the adjacent Riemann sheets, so they have disappeared from the physical illustration (c). 
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a < 0.5, then both ^i and (^3 are now on secondary sheets. This motion for variable a is 
shown as the dashed hne in the figure, and through the insets (a) to (c). 

The turning point, (t — Cij is real for all values of a and A. When A — 0, (i lies on top 
of the stagnation point, and as A — >■ 00, Ci moves leftwards along the real (^-axis towards 
the corner point by passing through |C| = cxd and traversing the free-surface. There are 
two special points: A = I (when (i has reached infinity upstream), and A = (b/a)^'^ (when 
Ci has reached infinity downstream) . This is illustrated for the case of a rectangular step 
in Figure m and insets (d) to (f). 



6.1. Inner analysis and reduction to the Airy Equation 

Turning points are important for two reasons: first, a turning point which lies on the free 
surface, as in the case oi A > 1, produces a change in one or both gravity and capillary 
waves from constant-amplitude oscillations to exponential decay [c.f. ^{dx/dw) in p.8p ]: 
second, turning points also generate Stokes lines and can thus lead to the birth of new 
exponentials. Recall from (I4.76P that, upon crossing the (g) > © or © > © Stokes line 
from C = C*) the base series switches on the exponential 



^cxp 



2^iQ±(C) 



exp 



1 /"^^ 



t< 



exp 



1 



^ dx± 



fdC 



(6.3) 



which we have written so that the integration passes through an arbitrary turning point, 
C = Ct. The Stokes lines from the turning points are given at locations where one expo- 
nential {e.g. gravity) reaches peak exponential dominance over the other (e.g. capillary). 
Thus, 



©>© 



©>© 



5R 



^ 



dx 



Ct dC 

^ dx 

Ct dC 



rdC 



ydC 



>5R 



^SR 



dx 



Ct dQ 
< dx 

Cr 



f < 



dC 



rdC 



(6.4) 



As ^ — > ^T, it is easy to verify from (16. 2p that the turning points are simple, and so we 
may write 

A^DiC-Cr), (6.5) 

where Z? e C Also, from (|3.14p and p.isp . we have that 

B± 



Qd 



(C-Ct)i/4 



(6.6) 



where B± is constant. We now write (|6.3p in a form that makes the connection to the 
Airy equation apparent. We introduce 



i^±(C)^-^(C-CT)^/*Q±exp 



1 



c* 



'^ ^ dC 

dC 



(6.7) 



so H± tends to a constant near the turning points. We also split the real and complex 
parts of x in (|3-10P with dx±/d( — {Si ± iS2), where 



Si{0^-{dw/dO'-^ and S2{0 = -{dw/dO^ 



(6.8) 



and dw/dC, = — 1/C is given by (|2.4p for the particular case of a channel geometry. The 



expression for qc 



New gravity- capillary waves (Part 2) 
in (16.31) can now be written as 



iexp 



H± cxp — 



C 



SidC 
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(6.9) 



Following the discussion of ^ we shall assume that ( = C,t '^s a first-order turning point 
(although this may not be true of all possible geometries specified by qo). If we like, 
we can re-scale q and C near C = Ct in order to remove the square-bracketed pre-factor 
of (16.91). The governing equa tions (|2.2ap - (j2.26p would then reduce to an Airy equation 
(c./. bender fc Orszaeill978l ). for which we can perform the local Stokes line analysis. 
However, it is faster to apply the following shortcut: using the local expression for A in 
(|6.5p . and combining with ()6.7p - (|6.9p gives 



gexp '^ h±{Ct) 



.T^iC-CTf^' 



(C - Ct)i/4 ' 

valid in the limit (^ — > Ct, where p e C is a constant that depends on ^t- 
We then note that the Airy switching requires 



(6.10) 



Turning Point 6^ 



iM- 



(C-Ct)i/4 ©^© ■ {C-CtV^^ ' ^^'^ (C - Ct)i/4 ' ^^-^^^ 

that is, the domi nant wave switches on the subdominant wave with a Stokes multiplier 
of i (| White! I2OO5I Chap. 11). The quantity M — ±1 inserts a sign that depends on the 
direction of analytic continuation past the Stokes line (the standard Airy equation, for 
example, has M = 1 for analytic continuation in the counter-clockwise direction, relative 
to the turning point). 

In terms of the outer variables, g, this shows that as we analytically continue across a 
Stokes line from the turning point at ^ = (t, the already existent capillary/gravity wave 
in (jO]) 

r<dx. 



^cxp 



27riQ-i 



X exp 



C* 



dC 



dC 



(6.12) 



switches on a gravity/capillary wave, which satisfies the transition rule: 



Turning Point . _ _ 

goxp — ^^ ^ — > 9cxp + «M 



© 



2TTiQ^ 



e'< 



X exp 



Ct 



dx± 
dC 



dC 



C 



B^ 

dx^ 



. # ''' 



X exp 



dx^ 



<.#''^' 



(6.13) 



where B± is given by (|6.6p . We have written the integrals in (|6.13p so as to emphasise 
the turning-point process. Notice that the last exponential in (j6.13p clearly indicates the 
presence of a gravity/capillary wave [ (I6.12p with a negated exponential]. The first two 
integrals of (|6.13p . however, show that the original capillary/gravity wave was generated 
from C* on the ± branch, went around the turning point C.x, and returned to C* along 
the =F branch. The values of most of the quantities {e.g. Q±, B±, and M) will not be 
necessary for the following discussion, but can be numerically computed if required. 



7. New solutions for the rectangular step {a = 1/2) 

In the spirit of King fc Blooii ( 19871) and IChapman fc Vanden-Broeckl ( 2006[ ) , we now 
study the simplest case for a non-trivial step: the rectangular geometry. In Part 1 of our 
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Figure 3. Positions ol the two turning points (|7.ip as a function oi A — ^t / P for the case of a 
rectangular step with a = 1 and 6 = 2. The Ti turning point a,t C, = C,i is shown thick, and the 
Tb turning point a.t C, = C,2 is shown thin. 



work, it was shown that for the case of a hnearised step of small height, two solutions 
are possible at small Froude and Bond numbers: one with capillary waves upstream and 
gravity waves downstream, and the other with localised solitary waves with decaying 
oscillations in the far field. Furthermore, this bifurcation arises because the capillary and 
gravity Stokes lines coalesce. 

However, for the nonlinear step, there is the added subtlety in the form of the turning 
points. As we know from ^ turning points not only change constant-amplitude waves 
to exponential decay, but also switch on secondary gravity or capillary waves. For the 
rectangular step, the two turning points are given by (|6.2p . simplified to 



Ci 



A 



1-Va 



and C2 



A 



1 



(7.1) 



VA ' 

and these are shown in Figure [3] for the step with a — \ and 6 = 2. Immediately, we 
recognise that there are at least three regions of interest: (i) both turning points on the 
solid boundary for Q < A < 1; (ii) one point on the surface and one on the solid boundary 
for 1 < ^ < (b/a)'^; and both points back on the solid boundary for A > {b/aY . 

In addition to these three regions, we must also note the positions of the Stokes lines 
relative to the turning points. It will be useful to introduce labels by 



C 
G 

Ti and Ti 
T2 



Capillary Stokes line from stagnation point, ( = —b 
Gravity Stokes line from corner point, ( = —a 
Capillary/gravity Stokes line from first turning point, C = Ci 
Capillary/gravity Stokes line from second turning point, C = C2 



The underline, Ti, is used to distinguish when the first turning point lies on bottom 
boundary and the overline, Ti, for when it lies on the free-surface. 

Next, examine the six subplots on the right of Figure |4j corresponding to six regions 
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labelled O and ©. These subplots indicate the qualitative arrangement of the turning 
points and Stokes lines for various values of A > 0. Solution-types are then classified 
by a sequence of labels, indicating the events as we traverse the free-surface from left to 
right (in the downstream direction). For example, a sequence like 

T\-T2-C-G 

in Region 1 corresponds to traversing the free-surface in the downstream direction, and 
crossing (i) the Ti Stokes line, (ii) the T2 Stokes line, (iii) the capillary Stokes line from 
the stagnation point, and (iv) the gravity Stokes line from the corner. 
Similarly, a sequence like 

C -Tl-G 

in Region 4 corresponds to (i) crossing the C-line, (ii) crossing the Ti point on the 
free-surface, and then (iii) crossing the G-line. 

The sequence of subfigures in the right of Figure 0] thus indicates that in the low- 
Froude, low-Bond limit, there are six possible solution types. Examine now the /3T-plane 
shown in the left frame of Figure ID The extremum cases of O and © correspond to the 
standard linearised solutions of Part 1 and confirms that the gray region of Figure [H 
which corresponds to 1 < A < (b/a)^, is further subdivided into @ to @. Thus, this 
figure shows that the usual dispersion 'line' which differentiates Rayleigh's two scenarios 
actually contains a series of solutions @ to @; the line's thickness and these new solutions 
are a manifestation of the nonlinearity (and largeness) of the step. As the step-height 
tends to zero, the line thickness tends to zero and we recover the standard linear picture. 

Remember that the physical sketches on the right subplots of Figure |3] are only illustra- 
tive projections of the Stokes lines from the upper half-(^-plane onto the physical domain 
(see Figure 2 of Part 1). The reader should not forget that in crossing a Stokes line along 
the free surface, a wave contribution due to the singularities in the lower half-^-plane 
also switch on; indeed the sum of these two contributions was used to derive (|4.8p . 

In the following discussion, it is sufficient to only consider the wave contributions from 
the upper half-^-plane. We shorten the notation for the arguments of the exponentials 
in (j4.7b[) and (|6.13p by expressing the integrals relative to either the stagnation point, 
( = —a, or the corner point, ( — —b. Thus 

Xeap^4£|,< and ^..av = 4£|t,c'. (7.2) 

We also denote the pre-factors for the exponentials by a calligraphic letter with a sub- 
script, if applicable, which corresponds to an interaction with a turning point. This leads 
to the three symbols: 



A 
C 
Q 



An arbitrary wave amplitude to be determined 

A capillary wave amplitude generated from the stagnation point 

A gravity wave amplitude generated from the corner point. 



Here are three examples: 

(i) ge""-- (ii) ^5Yt. e""^- (iii) {A + C)5^e^-- 

Example (i) denotes a gravity wave with amplitude Q switched on after crossing the Stokes 
line from the corner. Example (ii) is constructed in a two-step process: an (arbitrary) 
capillary wave with amplitude A encounters a turning point and switches on a gravity 
wave with amplitude Aj^, then encounters the second turning point and switches on a 
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Figure 4. The six regimes of interest for flow over a rectangular step, as shown in the /3r-plane 
(left) for the case of a = 1 and 6 = 2. The gray strip is bounded by A = 1 and A = (b/a) , 
which thickens for larger and larger steps. The arrangement of Stokes lines and turning points 
is shown on the right. Stagnation and corner points are circular nodes and turning points are 
square nodes. The arrangement in Region 6 is complicated and will be discussed in SI7.6I 
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capillary wave with amplitude Aq=^rp . Example (iii) is constructed by beginning with a 
capillary wave with combined pre-factors A and C (the former an arbitrary quantity, and 
the latter generated from the stagnation point), and then switching on a gravity wave 
due to the Ti turning point. Our notation and the associated expressions will be made 
clear in the many examples and figures to follow. 

7.1. Region 1 with Ti ■ T2 ■ C ■ G 

Consider the case oi A < 1. For most (larger) values of r, as we analytically continue 
across the free-surface, the sequence of events is Ti • r2 • C* ■ G, which involves the T2 Stokes 
line crossing the C Stokes line. In fact, for extremely small values of r, the lines may not 
actually cross, and so the sequence is Ti • C • T2 • G; however, both cases are equivalent 
because, as we shall see, the turning points do not play any role in determining the free 
surface. 

In order to derive the correct wave expressions, we traverse the free-surface, from left- 
to-right, beginning with a wave that satisfies the upstream radiation condition {i.e. a 
capillary wave). Upon reaching the end, we apply the downstream radiation condition, 
and this process provides us with a complete set of fully-determined waves. If we start 
with an arbitrary upstream capillary wave Ae ""^^ , then the sequence of events is 

©>© 



©>© 

-5^^^(yl + C)e^- 
®>© 

^^^^^(yl + C)e^-+ge^--. 
@>© 

There are no capillary waves downstream, bo A + C — Q, and the final solution has 

V- Stag. Corn. ^ v 

©>© © >© 

where the pre-factors, C and Q, are given by the jump condition (I4.76p . with 

C = -^^^ and g = _^^!^ (7.3) 

and Q± is computed from (|3.14p and (|3.15p . For the rectangular step, it can be shown 
by matching inner and outer solutions that 7c = and 7,3 = 6/5. In the end, the free- 
surface resembles the standard linearised solution with capillary waves upstream and 
gravity waves downstream, but now, with a wave-free region near the centre, similar to 
what occurs in Part 1 of our work. The result is sketched in Figure [Sf a). Note that the 
Stokes lines from the turning points remain inactive. 

In some problems, the intersection of Stokes lines can le ad to more subtle e ffects, such 



as the higher-order Stokes Phenomenon (see for example. iHowls et a/.l 120041 ). A simple 
consistency check is to analytically continue in a closed path around the intersection 
point, and to see if the solution returns to its original value. In this case, we see that the 
T2 Stokes line remains inactive, and so the intersection of the ® > © and © > © lines 
are not a concern. This will not be the case for the solution in Region 6. 

7.2. Region 2 with T\ ■ T2 ■ C ■ G 

As T is increased and /3 held steady, the Ti turning point moves upstream along the 
solid boundary and reaches —00 at A = 1. Solutions in Region 2 have A > 1 and the 
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(a) Region O with T\-Tt.-C -G 




(b) Region © with TY • T2 • C • G 




'////////////////////////////////// 
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(c) Region © with T\ ■ C ■ G 

Figure 5. Sketch of gravity-capiUary solutions in Regions O to ©. Waves sketched with the 
larger amplitudes are switched on by the stagnation point or corner; waves sketched with smaller 
amplitudes are switched on by turning points. Note that the gravity and capillary waves differ 
in both amplitude and wavelength, but we have not illustrated this difference. 



turning point T\ has now moved onto the free-surface and travels downstream as t is 
further increased. The turning point on the free-surface also implies that upstream waves 
will now decay, but downstream waves will remain constant. We start with an arbitrary 
upstream capillary wave Ae^^'^^ , and the sequence of events is: 



Ae 



X,. 



Ti 



©ic© 
T2 



^ Ae 



X,, 



„x. 



^W 



„Xc, 



,x. 



^^ (A + ^j,^j,Je— -r ^T, 



®>© 

Corn. 



T1T2 ' ^1'-' ' •^Ti'- 

> {A + ^jyt2 + C)e'^=''p + {Aj^ + g)e^'' 



®>© 

The last transition is known by (|4.7&|1 . so Q is again given by ()7.3p . Since there are no 
capillary waves downstream, and C is provided by the jump condition in (|4.76p . then 






T1T2 



flC 



(7.4) 
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We can associate the pre- factor A with the stagnation point, writing 



-4(0 = 



etc 



A{C), (7.5) 



for some unknown function, A{C,). The pre-factor, Aj^^ of the first gravity wave switched 
on by the Ti turning-point is given by the transition rule (j6.13p . and this allows us to 
relate A to A—. Care must be taken to express the gravitational wave quantities with 
respect to the corner, instead of the stagnation point [see the definition of Xgrav in (|7.2p ]. 
Next, the pre-factor, AtjTj, , of the second capillary wave switched on by the T2 turning 

point is also given by (|6.13|) . which relates A to A—j. . Lastly, (I7.4p and (|7.5p allows A 
to be solved entirely in terms of Q± and thus all quantities can be determined. 

Therefore, far upstream we have a decaying capillary wave, while far downstream we 
have a constant gravity wave. Also, switched on after encountering the Ti turning point is 
a (doubly) exponentially small gravity wave, which continues downstream. Finally, there 
is an even smaller capillary wave which has been turned on by the gravity wave across 
the T2-line. This is illustrated in Figure [Sljb), and in summary, 

©^© ^' ©>© ^'^' ^' 

@>© ^1 @>© ^ ^1 

In contrast to Region 1, the T2 line is now active, but the intersection of the two Stokes 
lines remains unproblematic; we can still analytically continue in a closed path around 
the intersection point and return to our original value. 

We note that the solutions in this region are actually contained in a rather small 
section of /3r-space (as shown in Figure |4]) . Once Ti has crossed the point that marks 
the intersection of the T2 line and the free-surface, a bifurcation occurs and the T2 line 
tends to w = —00 without intersecting the free-surface; this can be verified by numerical 
integration. 

7.3. Region 3 with Ti ■ C ■ G 

In Region 3, the only difference to the previous region is that the r2-line tends to it; = — oo 
without intersecting the free surface. Starting with an arbitrary upstream capillary wave, 
^g-f cap ^ ^jjg sequence of events is now: 

^e^cap — Zi — y Ae^^''^ + ytTj^e^s'"" 

®>© ^^ 

-^^[A + C)e^-^ + {A^ + g)e^^-^. 
® > © ' 

The downstream radiation condition requires C — —A. Like the previous regions, (j7.3p 
provides the expressions for C and Q, and (I6.13P relates A to Aj^. The final result is 

_4gXcap — Zi — y Ae^"''^ + AttB^'^"'^ 
©^© ^' 



^*"^-^ ^5^e^-" ^^^^^^ (A^ + g)e' 



®>© ^1 ©>© 

Thus the only difference between the solutions in Regions 2 and 3 is that in the latter. 
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the doubly-exponentially-small capillary wave has disappeared from the free surface. A 
sketch of the solution is shown in Figure El^c) . 

7.4. Region 4 with C ■ Ti ■ G 

For A still larger, Ti eventually crosses the intersection point of the C-line and the free- 
surface, and then the sequence is then C • Ti • G. If we start with an arbitrary upstream 
capillary wave Ae^'^'-^ , the sequence of events is 

©>© 



Ti 



©^© 



> {A + C)e^''--^^ + (yl + C)j^e'^'5"- 



-^ [A + C)e^- + ((^ + C)^ + C?)e^-v . 

® >© 
Again, we cannot have capillary waves downstream, so ^ + C = 0, and the coefficients, 
C and C/, are given by (|7.3p . The final result is 

^ V- Stag. Corn. ^ y 

®>© © >© 

The solution is sketched in Figure[6l[a), and we see that it consists of a decaying capillary 
wave upstream and a constant amplitude gravity wave downstream. 

7.5. Region 5 with C ■ G ■ Ti 

Eventually, Ti passes the point where the G-line intersects the free-surface, and so down- 
stream from this point, the gravity waves switch from constant-amplitude to exponen- 
tially decaying oscillations. For solutions in this region, it is easier to analytically continue 
from downstream to upstream, and negate each of the transitions [since the rules ()4.76| 
and (|6.13p are written for continuation in the downstream direction] . If we start with an 
arbitrary downstream gravity wave ^e''^^'"^ , the sequence follows 

_4e-^srav ^_Il yl^e-^^p + yle^s-- 

©ic© ^' 

©>© ^1 

^^{-A^-C)e^-^ + {A-g)e^^-^. 

We impose the requirement that there are no gravity waves downstream, so ,4 — ^ = 0. 
The values Bs and Cc are given in (|7.3p . and A and Ajr are related through (|6.13p . The 
final result is given by (from left to right) 

The solution is shown in Figure IHJb), where we see it consists of two decaying capil- 
lary waves upstream, and a gravity wave downstream. The gravity wave, which is first 
switched on by the Stokes line from the corner, decays until it reaches the turning point, 
and then switches to a constant amplitude wave downstream. 

7.6. Region 6 with G ■ G 

Once A = (b/aY^^, the Ti point has reached w = oo and for larger values of A, the 
turning point begins to move from right-to-left along the downstream solid boundary. 

There is one issue which needs to be addressed: in Figure Eljc), the ® > © line 
intersects the © > © line at a Stokes Crossing Point or SCP. If we analytically continue 
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-( GRAVITY )- 
-( CAPILLARY )- 



vwvwv wv 




®>© 



7777777777777777777777777777777777 



V777777777777777777777777777777777r/, 

(a) Region © with C - lY • G 




7///////////////////////////////// 



7777777777777777777777777777777777? 

(b) Region © with C ■ G • lY 




777777777779^7777777777^777777777 



(c) Region © with G • G 

Figure 6. Sketch of gravity-capillary solutions in Regions to ©. The solution in Region © 
should correspond to the second of Rayleigh's linearised profiles, but is complicated by the 
two crossing Stokes lines. In Appendix [XJ we demonstrate the existence of a third Stokes line 
(dashed), which originates from a secondary singularity (black circle). 



in a path that encircles the SCP, then we find that in order to avoid an inconsistency, 
the base solution must switch-off a capillary wave somewhere along the path; in other 
words, there must be a @ > © line, which goes through the SCP. However, neither the 
corner, nor the stagnation point, produces such a Stokes line. 

This issue is addressed in Appendix |X1 where we discuss three key ideas: (i) there is a 
second singularity [marked by the black circle in Figure IHlJc)]. This singularity does not 
appear in any of the early terms p.2ap - (l3.2dj) . but represents a singularity in the late 
terms o/ the late terms. It lies on the adjacent Riemann sheet from the originally defined 
corner point, C = ") S'Hd can be detected by setting x = 0, where from 



^l 



2t(?o 



d(/3. 



The contour F begins sX Q — a^ crosses the branch cut from the turning point, Ti, and 
continues along the secondary Riemann sheet of a/A- 

In the appendix, we also establish that (ii) the secondary singularity produces a © > © 
Stokes line that crosses the SCP [shown dashed in Figure[6ljc)]. Finally, (iii) the portion of 
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the ® > © line that joins the SCP to the singularity must be inactive; this switching-off 
of a Stokes line is due to the higher-order Stokes Phenomenon. These more complicated 
issues of secondary singularities, Stokes crossing points, and the higher-order Stokes 
Phenomeno n, have been s t udied by oth e rs in a variety of diff e rent s i tuations [see for 
examp l e Berk et all (fl983): IHowIs e^ all (12004): lOlde Daalhuisl (|2004 : IChapman eTal 
(|1999I) : [chapman fc Mortimeij (|2005[ ) : IChapman et all (|200l1 )]. 



The solution is then 

C + QtA e^"- -5^^^ C/T,e^-p Sc^^n^ ^or^ g^X,,. 

-^J ®>©-L ©>© ®>© 

where the values of C and Q are given by (|7.3p . The pre-factor Qt^ reminds us that 
the capillary wave is switched-on by taking the Q gravity wave and crossing the Stokes 
line from the the turning point; its value can be computed from (j6.13p . These doubly 
exponentially small capillary waves are switched off downstream upon crossing the Stokes 
line from the secondary singularity. 

Decaying waves on either side of the step are exactly what we expect from the linear 
analysis; however, we see that the result is not quite the same as in the linearised theory 
due to the presence of the doubly-small capillary waves upstream. The free surface is 
shown in Figure [6fc). However, w e note that, as for the case of pure gravity waves (see 



Chapman fc Vanden-Broeckl 120061 ) . we would expect a distinguished limit in the small 



Froude, small Bond, and small step height limits. Thus relating the results here with the 
work in Part 1 is not likely to be a trivial problem. 



8. Discussion 

Now having reached the end of our work on the gravity-capillary problem, we have one 
final query: do these new waves truly exist? 

Our theoretical predictions form the first attempt of exploring the previously unknown 
region of low-Froude and low-Bond numbers. Our results posit the existence of six differ- 
ent families of waves for flow over a step (O to © in Figured]), and reveals that the usual 
dispersion curve separating linear solutions O and © widens if we include the nonlinear 
nature of the geometry. While this structure is only valid at low Froude and low Bond 
numbers, it also sheds new light on the general complexity of gravity-capillary problem, 
which has been freely acknowledged in the past. 

Many open questions remain, some of which may ultimately hold the key to detecting 
these new waves (or discounting their existence) . Throughout this work, we have provided 
a clear analysis of the local issues of the gravity-capillary problem. For example, we have 
shown how the emergence of a Stokes line depends on the local angle of the obstruction, 
or how the switching-on of waves near turning points can be predicted via an Airy-like 
transition. Although we have not shown the tedious calculations, the crucial pre-factor, 
A, is another local issue, a nd can be derived by matching inner and outer solutions near 
the singularities feee lTrinhll2010l for more details). But what about the global properties 
of the Stokes lines? 

For example, can a Stokes line emerge from a singularity along a secondary Riemann 
sheet, only to later return to the primary sheet and intersect the free surface? Our analysis 
in Appendix \^ attempted to answer this question by drawing a comparison to a similar 
problem containing multiple singularities and crossing Stokes lines; we then used these 
results to propose the likely free-surface configuration for the complicated situation of 
Region 6 of Figure [Gjc). However, a more rigorous treatment of these issues remains an 
open problem. 
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Equally fascinating is the myriad of configurations that now seems possible for free- 
surface gravity-capillary flows over different geometries. A classic question is to inquire 
whether there are special geometries for which leading-order wave cancellation occurs and 
where, for example, the waves produced by one singularity cancels the waves produced 
by another. Although these issues regardi ng variable geometries have been studied in the 
context of pure gravity or capillary waves ( Chapraan fc Vanden-Broeckll2002t iTrinh etcd. 



l2011tlLustri et aLll2012tlTrinh fc ChaDmanl2013Ql ). a catalogue of the wave configurations 
for different geometries in the combined gravity-capillary case is a subject of future work. 
Finally, the most significant question is why have these waves not been previously 
detected? From the perspective of numerical simulations, one plausible difficulty concerns 
the so-called radiation problem: it is unclear how the upstream radiation condition can 
be accurately impleme nted in a numerical scheme. This d ifficulty is highlighted in, for 
example, the works of iForbed ( 1983[ ). IScullenI (|l998'), and iGrandison fc Vanden-Broeck 
(|2006 ). For gravity-capillary flow past a large obstruction, the resultant solutions are 
influenced by the chosen treatment of the far field, and it is possible that inherent errors 
in current numerical methods make it difficult to differentiate between the structure of 
our new gravity-capillary waves and the typical linearised solutions. Do these waves truly 
exist? Our theoretical results suggest that they do; the search for such configurations, 
either in nature or in the digital world, forms the basis of ongoing investigation. 



Appendix A. On the crossing of Stokes lines 

We address the issue of the crossing of Stokes lines in the gravity-capillary problem by 
drawing an analogy to the same phenomenon which occurs in a modified Airy equation: 



ey 



cX/n 



xy 



This equation was chosen because it not only possesses an (Airy) turning point at x = 0, 
but also an additional singularity at x — a. Writing y = ve^'"^ gives 



e v" + 2ev' + (1 — x)v = 



1 



(Al) 



where we also require w — > as |x| — > oo. We can now perform the standard asymptotic 
analysis of (jA 1[) . expanding v = J2 ^""^n, and giving the first two orders as 



vo 



Vl 



and in general, for n ^ 2, 



1 



(1 — x){x — a) 
-2 



(a — x)'^{x — a) (f — xY{x — a) 



2 ' 



(A 2) 
(A3) 



l-x 



•-'n-2 



2<-l 



(A 4) 



As n — >■ oo, we expect the singularities at a; = and x — a to produce factorial over 
power divergence of the late terms. Substituting the ansatz 

V{x)T{n + -i) 



[X(a^)]" 



+7 
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into (|A4p gives (x')^ ~ ^x' + (1 — a;) = 0, or solving yields 

X = / 1 ± s^/^ ds. 



(A 5) 



There are clearly singularities at x = 1 and x = a, and we shall focus on the exponentials 
generated by the latter. At a Stokes line emerging from the singularity at a; = a, the base 
series can switch-on exponentials of the form, e~^' "^ with x being one of 



Xi 



X2 



^3/2 



,3/2 



,3/2 



,3/2 



(A 6) 

(A 7) 



If we use ® to denote the base series, and ® and @ for the two exponentials, then the 
Stokes Phenomenon describes the process in which ® > ® or ® > @. However, there is 
also a turning point at a; = 0, for which the associated Stokes lines can cause interactions 
between the two exponentials. If ® > @, then this produces an exponential with 



X3 



2^3/2 



-a3/2, 



(Ai 



since the Airy transition simply requires switching the branch of x^''^ . Similarly, if @ > 
®, this produces an exponential with 



a; + 1^3/2 _^_^ 1^3/2^ (A 9) 

o o 

-1 + i are shown in the left frame of Figure [T] Notice 



X4 

The Stokes lines for the case a = 
that the ® > ® Stokes line intersects the ® > ® line at a Stokes Crossing Point (SCP). If 
we analytically continue in a circle around the SCP, we see that a third Stokes line, with 
® > ®, is needed in order to avoid an inconsistency. From (jA8|) . this seems to require 
a singularity at the point x — x^, where X3 = 0; for a = —1 + i, the singularity lies at 
X* « —1.423 — 0.509i. There is a similar singularity at a; « 0.249 + 0.883z, where X4 = 0, 
which we associate with a ® > ® line. The curiosity, however, is that these singularities 
do not appear anywhere in the base series, ®, and so they do not seem to be associated 
with any eventual divergence. 

To explore this in more detail, we define the Fourier transform as 



d{k) 



and take the transform of (|Aip . giving 



w(a;)e''''' dx, 



■ 2; 2^ 



2efcu — 



dv 
dk 



= V2^e"'''H{k), 



where H{k) is the Heaviside function, and we have assumed that 3(a) > 0. Solving this 
equation yields 



V = \/2tt exp 
or once inverted. 



V = I exp 

' — oo 



ek'^ + ik 



exp 



eu — iu + iau 



du, 



efc + ik ~ ikx 



exp 



dudk. 
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Figure 7. Stokes lines in the complex x-plane, shown with two different choices of branch cuts 
(left/right). White circles denote the singularities at a; = a and a; = 1, black circles are the 
secondary singularities, and a square denotes the turning point. The Stokes lines from a; = 1 are 
not drawn. Branch cuts are shown as a striped line, and the dashed lines are used to indicate 
points on a secondary Riemann sheet. Background contours are for jxsl (left) or |xi| (right), 
with darker regions indicating smaller values. 



We rescale k = s/e and i ~ w/e, giving 

(j){s;x) 



where we have defined 



exp 



0(w;a) 



exp 



4>{v;a) 



dv ds. 



(A 10) 



IV 



iv + aiv. 



In the limit e — > 0, the dominant contributions of the innermost integral in (jA 10[) come 
from the end points at t; = and u = s, or at the two saddle points at v — —i ± i^fa. 
Thus (jA10[) can be approximated as a sum of integrals of the form 



where i\) is one of 



V'(s) 



/ A(s)e'^('')/' 


ds. 


J-oo 




is(a — x) 




lb 9 

— — + s + «s — ISX 




is^ P 
— :r + * + IS ~ isx - 
o 


1 

3 



a± 



2a3/2 



(All) 

(A 12) 
(A 13) 

(A 14) 



Approximating the integral (lAllI) using (jA12p recovers the base series (jA2p - (IA3|l . The 
second and third expressions for ijj have saddle points at s = — i ± i^/x. Using the second 
expression (JA13|) . we see that the integral (|All| produces contributions e~'<^/^ with 



X 



-i+.±^.3/.. 



One of these comes from the singularity at a; = 1; the other comes from the same 
singularity, and then going around the turning point. 
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Finally, approximating the integral at the saddle points using the third expression for 
tp in (I A 14p . gives four possible expressions for x- Two of these are directly from the 
singularity at a; = a, and return xi f^nd X2 from (|A 6[) - (jA7[) . The other two, returning X3 
and Xi in (|A 8[) ^ (|A9p . are generated by the same singularity, but involve an integration 
contour which goes around the turning point. This verifies that the missing Stokes line 
should be there, though it does not shed much light on the mysterious singularity, a; = x*. 



Let us go back to the expansion. We can write the n*"^ term, «„, from (IA4[) as 



n+l 2n 



flfc, 



■' ^^ (a;-a)'=(l-a:)™' 

1 m=0 ^ ' ^ ' 



(A 15) 



with flfc^m given by 

afe,m = 2(A: - l)afe_i^,„_i - 2(m - 2)afc^„j_2 - (^ - 2)afe_2,m-i - (^ - 2) afc-2,m-i 
+ 2(fc - 1)(to - 2)afc_i^,„_2 - (m - 3)afc^,„_3 - (m - 3)^afe^„_3. 

We believe that if this recurrence relation is solved, and the method of steepest de- 
scents is used to approximate the terms (jA Va\ as ri —> oo, then the mysterious point, 
X = a:* , would appear as a singularity in the late terms o/ the late terms, Vn- The s e 'sec- 
ond ary' singular i ties an d SCPs have been encountered in lChapman fc Mortimeij ( 20051) 
and lHowls et al\ (|2004l ). 

Before moving to the gravity-capillary problem, we clarify a labelling issue: the pre- 
scription of X, by X3 = in (|A8P is correct only if the (g) > @ Stokes line can be 
continued from the SCP to the secondary singularity, x*, along the same branch of x'^'^ 
(as it is drawn in the left frame of Figure [7]). However, a different choice of the branch 
cut (as in the right frame) shows that xi = is the correct equation. To correct this 
ambiguity, we use an alternative notation. Let X(fc) correspond to a point on the fc"^ 
Riemann sheet (associated with the two branches of x"^'^); that is, X{u\ is mapped to 
g7rife^3/2^ with fc = 0, 1. From \Kh\ . we then write xi as 



Xi = 



1 



g7rJ*;(6)^3/2 



ds. 



(A 16) 



where A:(s) = or 1, depending on which sheet the integrand is currently on. This makes 
it very clear that the secondary singularity, x*, is still given by xi = Oi but x* may either 



lie on the sheet with a(o)7 or the sheet with a(i). Thus, X'i iii (|A8[) is simply 



X3 




g7rife(s)^3/2 



ds. 



except that x is restricted to the same sheet as an), since (|A8[) requires principal branches 
taken throughout. Writing x\ iii the form (jA 16|) serves to emphasize the fact that the 
contribution xz is found by beginning at ain\ and then crossing onto the secondary sheet. 

A.l. Crossing Stokes lines in the gravity- capillary problem 

We now address the issue of the crossing Stokes lines in Region 6 of iJ7.6l For x in 
(|3.10p . there are three types of branch points, but only the ones from the turning points 
are relevant. Remembering that the corner of the step is at C = a, and changing to 
integration in C, using (|2.4p . we denote, for the gravity wave. 



Xg 



1(1) 



2rgo 



ds, 



(A 17) 
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where k[s) = 1 on the branch of a^i) and fc(s) = on the branch of a(o) (associated 
with the capillary wave). In addition to the two turning points, there is a logarithmic 
branch point from the factor of s^^, and two square- root branch points from the two 
singularities of go- The former type adds a constant to x upon crossing the cut, while 
latter type switches go to —go upon crossing the cut. If we wanted to fully explore the 
different sheets, it would be better to write C, — CikiMMMM)^ ^^ keep track of how many 
times we have gone around each of the five branch points. However for our purpose, only 
the single turning point which causes the intersection is important. 

Remember: the issue is that the crossing of the (g) > © and © > © lines in FigurelHl^c) 
also requires a ® > © line in order to avoid an inconsistency. Like the example of the 
preceding section, this missing Stokes line comes from a secondary singularity, C = C*) 
not present in any of the early orders. To see this, we write (|A17p in the form 



Xg = 




Srgo 



- ds, (A 18) 

s 



and we further restrict ^ to the k — 0'^ Riemann sheet. Thus, (jA ISp is an integral from 
a(i) to a(o) via the branch cut from the turning point, and then proceeds from a(o) to (. 
The second integral is the usual capillary integral (hence the capillary wave). This new 
representation is used to make Figure [3 the most important feature of the figure is the 
missing singularity and its missing Stokes line, which indeed intersects the SCP. 

In addition to the appearance of the secondary singularity, intersecting Stokes lines 
are typically accomp anied by yet anothe r sub tlety: that of the Higher - Order Stokes Phe- 
nomenon [see again IHowIs et al\ ( 20041 ) and IChapman &: Mortimeii ( 20051 )]. This phe- 



nomenon specifies that at a SCP, Stokes lines may themselves switch-off. Indeed, this 
must be the case because if we analytically continue around the SCP in the inset of 
Figure [6l^c), we see that both portions of the ® > © line (on either side of the SCP) can 
not be active. It remains to determine whether the portion of the Stokes line connecting 
the singularity to the SCP is active, or whether it is rather the portion from the SCP to 
the free surface. 

To address this issue, we note in Figure |5] that if the dashed Stokes line is followed 
from the SCP to the secondary singularity into the lower half-plane, then it changes from 
a @ > © Stokes line to a ® > © Stokes line across the secondary singularity {(^ = (*). 
In fact, the portion of the axis where — & ^ C ^ Ci: corresponds to the Anti-Stokes lines 
(where ® = © = ©) from the corner {( = —b) and turning point {C, = Ci)- Across 
the secondary singularity, then, the Stokes line has changed character abruptly, and we 
argue by analogy to studies of the Higher-Order Stokes Phenomenon that this can not 
occur unless (* is a turning point (which it is not). Thus, the portion of the Stokes lines 
connected to the secondary singularity is inactive. Finally, we add that in Figure [H the 
dashed Stokes line is the only line associated with secondary Riemann sheets that we 
have chosen to show. There are other Stokes lines we have not shown (and in fact, other 
secondary singularities), but these are not relevant to the free surface waves. 
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Figure 8. The underlying contours and shading represent \xg\ from (|A 18|l . with dark regions 
indicating small values. Values of /3 = 1 and r = 1.5 were used, and a = 1, b — 2 for the corner 
and stagnation points (circles). Turning points are i^i « —0.310 and ("2 « —1.290 (squares), and 
the Stokes lines were computed using (|5.1|l and (|6.4|l . Note that with the exception of the Stokes 
line from (2, all the Stokes lines that begin in the upper-half plane will eventually intersect 
1^ G R^. Branch cuts are shown striped. The secondary singularity, located at (^' ~ —0.435, is 
shown as a black circle, and its Stokes line shown dashed (note that this Stokes line, as well as 
the singularity, lies on a different Riemann sheet than the other components of the figure). The 
inset confirms that if we analytically continue around the SCP, beginning with the base solution 
upstream, then a portion of the new Stokes line must switch off across the SCP (gray, dashed). 
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